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ABSTRACT 

We present a joint analysis of near-ultraviolet (NUV) data from the GALEX mission 
and (optical) colour profiles for a sample of 7 Brightest Cluster Galaxies (BCGs) in 
the Canadian Cluster Comparison Project. We find that every BCG which has a blue 
O ' rest-frame UV colour also shows a blue-core in its optical colour profile. Conversely, 

BCGs that lack blue cores and show monotonic colour gradients typical of old elliptical 
^ i galaxies, arc red in the UV. We interpret this as evidence that the NUV enhancement 

' in the blue BCGs is driven by recent star formation and not from old evolved stellar 

populations such as horizontal branch stars. Furthermore, the UV enhancement cannot 
' be from an AGN because the spatial extent of the blue cores is significantly larger 

^ ' than the possible contamination region due to a massive black hole. The recent star 

. formation in the blue BCGs typically has an age less than 200 Myrs and contributes 

mass fractions of less than a percent. Although the sample studied here is small, 
we demonstrate, for the first time, a one-to-one correspondence between blue cores 
in elliptical galaxies (in particular BCGs) and a NUV-enhancement observed using 
' GALEX. The combination of this one-to-one correspondence and the consistently 

, young age of recent star formation, coupled with additional correlations with the host 

cluster's X-ray properties, strongly suggests that the star formation is fueled by gas 
cooling out of the intracluster medium. In turn, this implies that any AGN heating of 
the intracluster medium in massive clusters only acts to reduce the magnitude of the 
cooling flow and that once this flow starts, it is nearly always active. Collectively, these 
. , results suggest that AGN feedback in present-day BCGs, while important, cannot be 

^ . as efficient as suggested by the recent theoretical model by proposed by De Lucia et 

al. (2006). 

Key v^fords: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - 
galaxies: evolution - cooling flows - X-rays: galaxies: clusters 
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1 INTRODUCTION 

Elliptical galaxies provide a critical test bench for many sce- 
narios of galaxy formation because they follow tight rela- 
tions in their photometric (optical), chemical and dynami- 
cal properties. Within massive groups and clusters, elliptical 
galaxies populate a tight region of the colour-magnitude re- 
lation (CMR) (Bower et al. 1992) known as the red sequence, 
which is mainly driven by an underlying 'mass-metallicity' 
relation (e.g. CaroUo et al. 1993). The lack of significant 
change in the slope and scatter of the CMR (e.g. Stanford 
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et al., 1998), the slow evolution of colours (e.g. Saglia et al. 
2000; Ellis et al., 1997), and line strength indices of cluster 
early- type galaxies out to z ~ 1 (e.g. Bernardi et al., 2003) 
indicate that the bulk of their stars that comprise these sys- 
tems formed over a relatively short period of time at high 
redshifts {z > 3) (see Renzini, 2006 for a recent review), 
though the issue of where these stars formed and how they 
were assembled into the observed elliptical galaxies has yet 
to be conclusively resolved. 

The simplest model for the formation of elliptical galax- 
ies — especially the most massive ones — is the 'mono- 
lithic collapse' model (see e.g. Larson 1974, Matteucci & 
Tornambe 1987, Pipino & Matteucci 2004, MerUn & Chiosi, 
2006). In this model, the elliptical galaxies are believed to 
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have formed from homogeneous coUapse of intergalactic ma- 
terial, foUowed by a rapid, massive starburst, during which 
nearly all of the available gas is converted into stars that 
comprise the observed stellar content of the elliptical galax- 
ies today. In its simplest form, the monolithic model predicts 
that the vast majority of the present-day elliptical systems 
were already be in place by z ~ 3 and that once formed, have 
been evolving passively towards the present. Consequently, 
this model predicts that ellipticals galaxies should be largely 
red and dead today. 

Detailed recent observations of elliptical galaxies, how- 
ever, indicate a more complex formation history: Taylor et 
al. (2008), for example, find that at the most only 20% of the 
local red sequence galaxies with stellar mass Af* > 10^^ Mq 
were already in place by jz ~ 2 (but see Whiley et al. 2008, 
who do not detect any significant change in the stellar mass 
of the BCG since z ~ 1). Kormendy et al. (2008) argue that 
the mass-dependent variations in the structural properties 
of elliptical galaxies is inconsistent with the idea of a uni- 
form, synchronized formation history. And using rest-frame 
UV observation^ to probe the finer details of the star for- 
mation histories in elliptical galaxies, Kaviraj et al. (2007a; 
2008) find compelling evidence that while the bulk of the 
stars that comprise the elliptical galaxies are old, the galax- 
ies continue to form stars at a reduced rate over the lifetime 
of the Universe. Specifically, using GALEX (UV) and SDSS 
(optical) photometry, Kaviraj et al. (2007a) find that at least 
30% of nearby {Q < z < 0.11), massive (Mr < -21) ellipti- 
cals show unambiguous signatures of star formation within 
the last Gyr, contributing up to a few percent of the stellar 
mass of the galaxy. 

An alternative model for the formation and evolution 
of galaxies is the hierarchical clustering model in which the 
these systems are assembled via successive mergers of lower 
mass systems, with both the formation of the lower mass 
progenitors occurring at an earlier epoch and the merger 
rate being larger in regions of the Universe characterized 
by higher than average total mass density. The hierarchical 
model is an inescapable prediction of the currently favored 
cold dark matter models for the formation and evolution 
of large-scale structure in the Universe. Strictly speaking, 
the hierarchical model only speaks to how the mass is as- 
sembled, not to the details of star formation, which is what 
the observations are sensitive to. The latter depends on a 
series of complex, highly nonlinear baryonic physics, includ- 
ing the manner in which galaxies acquire their gas (Keres 
et al. 2005), the local environment where the galaxy and 
its progenitors reside, the details of the cooling and heating 
processes that the gas is subjected to (e.g. Silk 1977; White 
& Rees, 1978; Scannapieco et al. 2005; Kaviraj et al. 2005, 
2007b), etc. The most recent detailed study of the forma- 
tion history of elliptical galaxies by De Lucia et al. (2006) 
shows that while the bulk of the stars that make up the most 
massive ellipticals form rapidly and high redshifts, they do 
so in a number of different progenitor systems and typi- 
cally, the assembly of the progenitors into a single massive 



^ In contrast to the optical spectral ranges, the UV is highly 
sensitive to even small fractions of young stars (younger than 
about a Gyr) , making it an excellent probe of the low-level recent 
star formation. 



elliptical galaxy occurs much later. Whether the observed 
structural, chemical and photometric properties of massive 
elliptical galaxies is compatible with hierarchical late-time 
assembly remains to be seen (see, for example, Cimatti et 
al. 2006; Ciotti et al. 2007; Pipino & Matteucci 2008; Pipino 
et al. 2009, Maiohno et al. 2008, and Whiley et al. 2008 for 
observational arguments against late-time assembly). 

In the context of the hierarchical models, the most im- 
portant of the recent innovations is the inclusion of AGN 
feedback. AGN feedback plays a central role in driving 
gas out of the progenitor systems, thereby both truncat- 
ing the initial starburst in these systems and ensuring that 
any subsequent mergers involving these objects are gas-poor 
(i.e. "dry mergers"). AGN feedback is also critical for pre- 
venting the cooling of any gas that subsequently accumulates 
in the halos of the massive elliptical galaxies or the larger 
systems (e.g. groups and clusters) that they are embedded 
in. Prior to the inclusion of AGN feedback, the hierarchi- 
cal model produced elliptical galaxies that were both much 
too luminous and much too blue compared to the present- 
day population of massive elliptical galaxies. Unfortunately, 
details of how AGN feedback actually works is not well un- 
derstood and the properties of the elliptical galaxies in the 
hierarchical model depends sensitively on how AGN feed- 
back is modeled. 

In this paper, we will use recent observations of a spe- 
cial class of massive elliptical galaxies, the Brightest Cluster 
Galaxies (BCGs), to gain additional insights into processes 
shaping the formation of elliptical galaxies. We are especially 
interested knowing whether these systems show evidence of 
late-time star formation, the conditions that facilitate late- 
time star formation and in understanding the nature of AGN 
feedback in these systems, at least at the present-time. As 
the most massive, most luminous galaxies in their local en- 
vironment these galaxies have long been of special inter- 
est. In spite of their special status, they have much in com- 
mon with other massive early-type galaxies. For example, 
the mean stellar ages and metallicities of BCG galaxies are 
similar to those of non-BCG ellipticals of the same mass 
(Fisher et al., 2005; Brough et al. 2007, von der Linden et al., 
2007). More importantly, within the context of the hierar- 
chical model, the essential features of the formation history 
of the BCGs is predicted to be similar to those of massive 
non-BCG early-type systems. Specifically, while nearly 80% 
of the stars that comprise the BCGs form at z > 3, they 
form in several distinct systems. The BCGs themselves do 
not take on a distinct identity until after z ~ 0.7. (De Lucia 
& Blaizot 2007; Romeo et al. 2008). 

The BCGs that emerge from the theoretical modeling 
by De Lucia & Blaizot (2007) are, for all intensive pur- 
poses, red and dead. They experience virtually no star for- 
mation after the initial burst because the modeling assumes 
an extremely efficient form of AGN feedback. However, sev- 
eral recent studies (non-UV based) have reported examples 
of ongoing star formation in the brightest cluster galaxies 
(Cardiel et al. 1998, Crawford et al. 1999, Edge 2001, Goto 
2005, McNamara et al. 2006, Wilman et al. 2006, O'Dea et 



It is often incorrectly assumed that BCGs always reside at the 
center of the cluster potential. As shown by Bildfell et al. (2008), 
this is not necessarily so. 
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al. 2008, Bildfell et al. 2008, Cavagnolo et al. 2008, RafTerty 
et al. 2008); most of these BCGs reside in cool core clusters. 
Additionally, Hicks & Mushotzky (2005) noted an excess in 
the UV flux — as determined from the XMM-Newton Opti- 
cal Monitor — in many (but not all) the cooling flow clusters 
in their sample, which they interpreted as evidence for star 
formation. Most recently, Bildfell et al. (2008) undertook a 
comprehensive study of 48 clusters that span a wide range of 
X-ray characteristics. Specifically, they analyzed the surface 
brightness and color profiles of the BCGs hosted by these 
clusters, seeking to relate the resulting trends to the relative 
location of the BCGs within the clusters as well as to their 
global X-ray luminosity (Lx) and temperature (T^). They 
found that 25% of their BCGs had bluer colors in their cen- 
tral regions, which they interpreted as evidence for ongoing 
star formation. They also found that these blue core systems 
only occurred in cool core clusters (see also Cavagnolo et al 
2008) and then only if the BCG is located at the cluster 
center. For completeness, we note that cool core clusters are 
systems whose central gas density is sufficiently high that 
the corresponding radiative cooling timescale is < 5 Gyrs. 
(See McCarthy et al. 2004; 2008 for a detailed discussion of 
the diversity in the cluster population, including the distinc- 
tion between cool core versus non-cool core clusters). 



The mounting evidence for active star formation poses 
a challenge for models that invoke strong AGN feedback. 
Bildfell et al (2008) suggest that in the systems that they 
studied, heating by AGN feedback may be offsetting most of 
the radiative losses suffered by the hot gas surrounding the 
BCGs but not all. Therefore, the gas cools but at a signifi- 
cantly reduced rate. If this is indeed the case, it represents an 
important clue into how AGN feedback operates. Our long- 
term goal, therefore, is to understand that implications of 
ongoing star formation in BCGs within the broader context 
of galaxy formation. As a first step, we need to establish that 
the blue core phenomenon is indeed associated with star for- 
mation. As discussed below, we do so by demonstrating that 
the blue core are unambiguously linked to UV-enhancement. 
Next we quantify the extent of recent star formation. Specif- 
ically, we are interested in determining the age of the last 
star formation event and the fraction of the total stellar 
mass that it gave rise to. Finally, we review the conditions 
characterizing the environments of the BCGs in our sample 
and use the apparent correlations to identify the source(s) 
of the cold gas that is fueling the residual star formation. 



In order to carry out our analysis, we start by cross- 
matching the 48 BCGs in the Bildfell et al. catalog with 
archival data from the GALEX GR3. We find 10 BCGs in 
common, of which we analyze 7. In Sec. [2l we briefly sum- 
marize the main characteristics of the Bildfell et al. sample 
and the data retrieved from the GALEX archive. In Sec. |3l 
we present the model set up used to estimate the character- 
istics of the recent star formation (ages and mass fraction) 
in each BCG. Finally, our results are discussed in Sec. ID 
and the conclusions are presented in Sec. [3 The cosmolog- 
ical parameters used throughout this work are Ho ~ 70 km 
s"^ Mpc"\ = 0.3 and = 0.7. 



2 DATA 

The galaxies studied by Bildfell et al. (2008) are drawn 
from the Canadian Cluster Comparison Project (CCCP) 
The CCCP is a study of an X-ray selected sample of 50 
clusters from the Advanced Satellite for Cosmology and As- 
trophysics (ASCA) catalog of Horner (2001), with redshifts 
in the range 0.15 2 ^ 0.55 and most clusters having X- 
ray temperatures Tx > 5 keV. Almost all clusters have also 
either Chandra or XMM-Newton high resolution data. The 
CCCP clusters span the observed range in scatter of clus- 
ter X-ray and S-Z scaling relations, and thus create maxi- 
mal leverage for the investigation of baryonic feedback ef- 
fects. The combination of X-ray and deep, multi-filter op- 
tical data makes the CCCP an excellent sample for study- 
ing the links between BCGs and their host cluster environ- 
ments. The deep optical imaging data was obtained using 
the Canada- France-Hawaii Telescope (CFHT). The full sam- 
ple comprises a set of 30 clusters observed in g' and r' filters 
using the MegaCam detector and 20 clusters observed in the 
B and R filters using the CFH12K detector which were taken 
from the CFHT data archive at the Canadian Astronomical 
Data Centre (CADC). 

A quarter (13/53) of the BCGs in the CCCP show 
colour profiles that become increasingly blue in the core. 
These galaxies are displaced from their host-cluster optical 
red sequence by ~ 0.5 to 1.0 mag in {g'-r'), which is im- 
portant for optical cluster studies that may reject the BCGs 
based on colour. The blue cores in these systems cannot 
be explained by AGN point-sources because (a) they are at 
least twice the size of the seeing disk and (b) Type II AGN 
do not contaminate the optical spectrum at more than a few 
percent level (Zakamska et al. 2006a; 2006b). The above ar- 
guments also apply to the NUV, in that the contamination 
due to AGN is less than 15% in the UV (see Salim et al. 
2007 for a more detailed discussion). Moreover, Kauffmann 
et al. (2007) examined a sample of nearby AGN hosts and 
find that their UV emission too is quite extended, thus un- 
likely to emanate from the central AGN, whereas it seems to 
be associated with star forming discs. For these reasons, we 
shall assume that all of the NUV fiux is purely stellar in ori- 
gin and compute star formation rates accordingly. Strictly 
speaking, however, the star formation rates that we derive 
represent upper limits. 

The Bildfell et al. sample was cross-matched with pub- 
licly available UV photometry from the second data release 
of the GALEX mission (Martin et al. 2005). GALEX pro- 
vides two UV filters: the far-ultraviolet {FUV), centered 
at ~ I530A and the near- ultraviolet (NUV), centered at 
~ 23IOA. Note that, since the bulk of our sample are at in- 
termediate redshifts {z > 0.15) we only use the NUV filter in 
this study, since the NUV filter traces the spectrum between 
rest-frame NUV and FUV. The cross-matching produced 10 
BCGs which have at least a detection in the NUV filter. 
The remainder of the galaxies are not observed partly be- 
cause of a lack of sufficient spatial coverage in Galex data 
and, possibly, short exposure times in some Galex fields. 
The positional matching is performed within the GALEX 
fiducial angular resolution of 6". NUV magnitudes and ex- 
posure times for the GALEX fields are given in Table (TJ 
along with other information related to the galaxies such as 
redshifts, J2000 coordinates, effective radii, and the type of 
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Table 1. Galaxies from Bildfell et al. (2008) present in GALEX DR2 



BCG Name 




otBCG 


<5bcg 


r',R g',B r^' 


Core 


R* 


NUV (2300A) 


exp 


rcf 




(J2000) 


{J2000) 


mag mag kpc 


colour 


kpc 


mag 


ks 



Abell 1835 
Abell 1942 
Abell 2111 
CL0910+41 
MS0440+02 
MS0451-03 



0.25 
0.22 
0.23 
0.44 
0.19 
0.54 



14 01 02.10 

14 38 21.88 

15 39 40.52 
09 13 45.52 
04 43 09.92 
04 54 10.84 



+02 52 42.69 
+03 40 13.34 
+34 25 27.46 
+40 56 28.54 
+02 10 19.33 
-03 00 51.39 



15.97 
16.35 
17.26 
18.93 
15.06 
19.38 



17.17 
17.38 
18.50 
18.81 



20.75 45.5 



+0.2 
0.2 

26.4t°:l 

+0.5 
0.4 

102.0+Oj 

+0.3 



43.7"* 



32.6"* 



-0.6 



blue 

red 

red 

blue 

red 

blue 



19 ±2 



40 ±8 



20 ±2 



19.13 ± 0.02 
22.68 ± 0.28 
23.41 ± 0.35 
19.58 ± 0.11 
21.84 ± 0.44 
20.65 ± 0.28 



1.7 
1.7 
7.3 
0.4 
0.09 
0.11 



7 


Abell 2219 


0.2256 


16 40 19.85 


+46 42 41.30 


15.98 


18.40 


61 7+°-'' 


red 




22.66 ± 0.27 


2.75 


8 


Abell 2390 


0.2280 


21 53 36.84 


+17 41 44.10 


17.33 


19.22 


18 4+°-^ 


blue 


25 ±5 


20.94 ± 0.07 


1.96 


9* 


CL0024+16 


0.39 


00 26 35.68 


+17 09 43.48 


17.45 




60 1 +^-^ 


red 




21.91 ± 0.09 


15.42 


10 


MS1512+36 


0.3727 


15 14 22.51 


+36 36 21.30 


17.66 


19.94 


57 6+° '' 
"'•"-0.9 


blue 


15 ±2 


21.40 ± 0.34 


0.07 



COOP targets from Bildfell et al. (2008). The first group of six BCGs was observed with CFHT MegaCam, while the second group 
was observed with the CFH12K camera. The galaxies marked with a * have not been used for the analysis (see Section 2). Redshifts 
(z) listed are spectroscopic. The (r', g') and {R, B) apparent magnitudes are in AB and Vega system, respectively (see text). The 
associated errors are typically 0.03 mag. Note that for the reasons outlined in section [2] no g' or B magnitudes were obtained for 
CL0024+16 and MS0440+02. The r' best-fit surface brightness profile parameter is given to be compared to the size of the blue 
region ij» - if a blue core is present. Core colour is given as 'red' or 'blue' based on the shape of the inner colour profile (see Bildfell et 
al., 2008). Galex NUV magnitudes are in AB system. Galex NUV exposure times (exp) in units of 1000s. 



core (red or blue) as described in Bildfell et al. In each case 
the GALEX image with the longest exposure time was used. 

We note that, in general, galaxies harbouring recent star 
formation are more likely to be detected than non recent star 
formation-galaxies (since their UV flux will be higher), es- 
pecially at higher redshift. Thus at any redshift the chances 
of detecting a UV blue BCG is higher than its UV red coun- 
terparts. However, since the goal of this work is to verify 
if a strong NUV flux is unambiguously linked to blue-core 
galax;ies (and vice versa) the incompleteness in the detection 
of UV red BCGs does not affect our results. 

Out of the 10 BCGs that are detected by GALEX, 5 
have blue cores and 5 have normal (red) cores. However, this 
sample reduces to 7 objects because we have to exclude the 
BCGs in MS0440, CL0024 and MS0451 for the reasons out- 
lined below. The MS0440 and CL0024 BCGs are excluded 
because there are systematics due to crowding in both of 
these systems. Moreover, the MS0440 BCG lies right at the 
edge of a very bright stellar reflection artifact, which causes 
problems for proper background subtraction. 

In the case of the MS0451, there is a large foreground 
spiral that overlaps with its BCG. To minimize contami- 
nation, Bildfell et al. were required to mask a significant 
fraction of the BCG before measuring a total magnitude 
from its optical image. The central regions of the BCG are 
more heavily masked than the outer regions, which leads the 
colour to be dominated by that of the outer regions making 
the galaxy optically red. Since we extract the UV photom- 
etry directly from the GALEX GR3 database (without any 
masking) and given the large PSF of GALEX (~ 6 arc sec- 
onds) it is likely that the star formation in the overlapping 
spiral contributes to the UV signal we measure. This results 
in the MS0451 BCG having a very red {g' — r') colour and 
an extremely blue UV-optical colour, rendering it difficult to 
perform accurate parameter estimation using realistic mod- 
els. 



3 MODELLING THE SPECTRUM: 

QUANTIFYING THE AGE AND MASS 
FRACTION IN YOUNG STARS 

We estimate parameters governing the star formation his- 
tory (SFH) of each galaxy by comparing its GALEX NUV 
and CFHT optical photometry to a library of synthetic pho- 
tometry, generated using a large collection of model star 
formation histories. Speciflcally, we follow Kaviraj et al. 
(2007b) and make use of g'-r' & (NUV-r) colours, with the 
near-UV constraining the amount of recent star formation 
and the red optical light providing the normalization. 

Each model SFH is constructed by assuming that a 
burst of star formation at high redshift {z ~ 3) is followed by 
a second burst, which is allowed to vary in age between 0.001 
Gyrs and the look-back time corresponding to z = 3 in the 
rest-frame of the galaxy, and mass fraction between and 
1. Both bursts are assumed to be instantaneous. A Salpeter 
(1955) initial mass function is adopted for each burst. Figure 
[T] shows a schematic representation of the model SFHs. As 
an alternative, we could have modeled the second star for- 
mation event as extended in time (or even, continuous). We 
are unable to do so, however, because this approach requires 
more wavelength coverage than currently available in order 
to break the degeneracies between age, metallicity and dust 
content. 

To build the library of synthetic photometry, each 
model SFH is combined with a single metallicity in the range 
O.IZ0 to 2.5Zq and a value of dust extinction parametrised 
by Eb-v in the range to 0.5. Photometric predictions 
are generated by combining each model SFH with the cho- 
sen metallicity and Eb~v values and convolving with the 
stellar models of Yi (2003) through the GALEX NUV and 
CFHT optical flffers. The model library contains ~ 750, 000 
individual models. Finally, since our galaxy sample spans a 
range of redshifts, equivalent libraries are constructed at the 
redshift of every galaxy. 

The primary free parameters in this analysis are the age 
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Figure 1. Model SFHs (see Section 3) are constructed by assum- 
ing that an instantaneous burst of star formation at high redshift 
(z = 3) is followed by a second instantaneous burst which is al- 
lowed to vary in age and mass fraction. The main free parameters 
are the age (tyc), mass fraction (fyc) of the second burst, tyc is 
allowed to vary from 0.001 Gyrs to the look-back time correspond- 



ing to z - 
and 1. 



: 3 in the rest-frame of each galaxy, fyc varies between 



0.10 Gyrs + 
0.50 Gyrs SK 
1.00 Gyrs O 
2.00 Gyrs A 



Figure 2. The sensitivity of the UV to young stars. We assume 
two instantaneous bursts of star formation, where the first burst 
is fixed at z = 3 and the second burst is allowed to vary in age and 
mass fraction. The near-UV (NUV) colour of the composite stellar 
population is plotted as a function of the age (symbol type) and 
mass fraction (x-axis) of the second burst. It is apparent that 
even a small mass fraction (~1%) of young stars 0.1 Gyrs 
old) causes a dramatic change in the NUV — r colour compared 
to what might be expected from a purely old stellar population 
{NUV - r ~ 6.8). 



(tyc) and mass fraction (fyc) of the second burst, which pro- 
vides the young component ('yc') of the stellar population 
in the galaxy. For each galaxy, parameters are estimated by 
comparing each observed galaxy to every model in the syn- 
thetic library, with the likelihood of each model (exp — 
calculated using the value of computed in the standard 
way. From the joint probability distribution, we marginalise 
over the dust and metallicity to extract the tyc vs. fyc prob- 
ability space (see next section) for each galaxy. This space 
shows, for every galaxy, the probability that the recent star 



formation is described by a particular value of tyc and fyc- 
This allows us to explore the characteristics of the recent star 
formation in the individual BCGs and correlate them to the 
presence (or absence) of a blue core in the optical image. The 
error in the is calculated by adding, in quadrature, the 
observational uncertainties and typical errors in the stellar 
models, which we assume to be 0.05 mag in each optical filter 
and 0.1 mag for the GALEX UV fiher (Yi 2003). Note that 
the method is similar in design to past techniques designed 
to detect small 'frostings' of young stars over a underlying 
old stellar population, either using rest-frame UV data (e.g. 
Ferreras and Silk 2000) or spectroscopic line indices (e.g. 
Trager et al. 2000). 

The sensitivity of the UV to young stars is demon- 
strated in Figure [S] We assume two instantaneous bursts 
of star formation, where the first burst is fixed at z = 3 and 
the second burst is allowed to vary in age and mass fraction. 
The near-UV colour of the composite stellar population is 
plotted as a function of the age (symbol type) and mass frac- 
tion (x-axis) of the second burst. It is apparent that even a 
very small mass fraction (~1%) of young stars (~ 0.1 Gyrs 
old) causes a dramatic change in the NUV — r colour com- 
pared to what might be expected from a purely old stellar 
population {NUV — r ~ 6.8) or from the analysis of any 
optical colour-magnitude relation. Given that typical obser- 
vational uncertainties in the NUV — r colours from modern 
instrumentation are ~ 0.2 mag, the usefulness of the UV in 
detecting residual amounts of recent star formation becomes 
quite apparent. Note that the NUV filter is taken from the 
GALEX filter set. 



4 RESULTS AND DISCUSSION 

In Figure|31 we present the UV-optical colours of the 7 BCGs 
in our final sample. In Figures |4] and [5] we present the quan- 
titative SFH parameter estimation in terms of t^ and fyc 
for the red-core and blue-core BCGs respectively|j 

Since the BCGs span a large range in redshift we do 
not simply show their {NUV — r) colours because their 
'red sequence' positions (i.e. the position of a old, passively 
evolving population) are very different. Instead, we show, in 
Fig- El the difference between the {NUV — r) colour of each 
galaxy and the 'red sequence' position at its given redshift, 
calculated using a dustless, solar-metallicity simple stellar 
population (SSP) forming at 2 = 3. The local early-type 
population (Kaviraj et al. 2007a) is shown over-plotted us- 
ing the small black dots. The value of Mr for the BCGs pre- 
sented here is calculated by converting from the magnitudes 
in the CFHT Mould and Megaprime r-band filters (from 
Bildfell et al.) using a solar-metallicity SSP that forms at 
z = 3. Note that this local population is drawn mainly from 
the field and its local density distribution varies strongly 
from that of the BCG population (see Figure 5 in Schawin- 
ski et al. 2007). 

Two results are immediately apparent on inspection of 
Figs. [SHS] First, the BCGs that do not exhibit blue cores 



We caution the reader that the estimates of tyc and fyc for the 
BCG in MS1512 are not particularly well constrained because of 
the large error in NUV flux due to extremely short exposure — 
see Table 1. 



6 Pipino et al. 




.000 



o 0.100 



0.010 



0.001 




10" 



10" 



10° 



Figure 3. The difference between the red sequence position of 
each BCG at its given redshift, calculated using a solar metal- 
licity simple stellar population (SSP) which forms at 2 = 3, and 
the {NUV — r) colour. The local early-type population is shown 
over-plotted using the small lack dots. The numbers identify the 
galaxies according to the order in Table 1. 



(i.e. nos. 2, 3, 7) are much closer to the NUV color of their 
respective red sequence positions: d{NUV — r) < 1.5. In 
other words, the red optical color corresponds to red NUV 
color. This is not necessarily surprising but it does provide 
us with a baseline. What is interesting, however, is that all 
of these UV-red BCGs are somewhat displaced by from their 
respective red sequence positions, suggesting that the star 
formation history of even these UV-red BCGs is inconsistent 
with their stellar content having formed entirely at high red- 
shift {z 3) and that they likely experienced low-level star 
formation events within the past ~ 2 Gyrs. This is the 'time 
window' during which the UV signal remains detectable (see 
top panel of Fig. 7 in Kaviraj et al. 2007a). This conjecture 
is borne out by the fact that the {tyc, fyc) map for the red- 
core BCGs (Figure|4ll is not concentrated exclusively at very 
old ages. In the case of A2111 BCG, for example, the recent 
star formation is consistent with intermediate age popula- 
tions (since tyc is around 1 Gyr and fyc is around a few 
percent). 

Given our small sample, we cannot reject the possibility 
that the d{NUV — r) for our red core BCGs is likely biased 
high due to the requirement that they must be detected in 
the UV before entering the sample. Based on this, we would 
expect that any missing red core systems would likely lie 
even closer to the red sequence than the ones shown here, 
further emphasizing the gap in UV color between the opti- 
cally blue core BCGs and their red core counterparts. On 
the other hand, the offsets exhibited by our red-core BCGs 
appears to be comparable to those of the local early-type 
population as a whole and suggest that low-level star forma- 
tion events are not uncommon or restricted to BCGs only. A 
careful analysis of the UV light in a large sample of elliptical 
galaxies indicates that luminous elliptical galaxies form up 
to 10-15 % of their stellar mass after z ~ 1, although for 
the bulk of the population, the typical young stellar mass 
fraction is much smaller, around a few percent (Kaviraj et 
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Figure 4. Probability maps for tyc and fyc for the UV-red BCGs. 
Recall that tyc and fyc represent the age and mass fraction of 
the recent star formation in the galaxy. From top to bottom the 
BCGs shown are: A1942, A2111 and A2219. Wc find that the 
'young component' ('yc') in each of these galaxies is composed of 
intermediate age populations with ages greater than 1 Gyr. Note 
that none of the galaxies are consistent with purely old (> 3 Gyr 
old) populations. This is consistent with the fact that all the UV- 
red BCGs lie at least 1 mag away from the red sequence positions 
at their given redshifts. 
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al. 2007a; 2008; see also Kaviraj 2008 for a review of these 
recent results). 

Second and perhaps most interestingly, Figs. [SHU re- 
veal an unambiguous one-to-one correspondence between UV 
'blueness' and a blue core m the optical image in our BCG 
sample. Every BCG that is far from the red sequence po- 
sition (i.e. has d{NUV — r) > 2.5), which is rarely seen in 
the regular ellipticals, is also a blue core system, and vice 
versa. Given that we associate the excess UV flux with star 
formation, this one-to-one correspondence, combined with 
the fact that UV-red BCGs do not show anomalous colour 
profiles, bolsters our earlier assertions (as well as those of 
Bildfell et al. 2008) that the blue-core BCGs show the fea- 
tures that they do because they host active star formation 
in their central regions. This association is further strength- 
ened by independent confirmation of ongoing star formation 
based on optical and infra-red spectroscopic work in three 
of our four blue-core systems. These three galaxies and the 
measured rate at which they are forming stars in their cen- 
tral regions are as follows: A2390 is estimated to be forming 
stars at ~ 5 Mq yr"\ CL0910 at ~ 40 Mq yr"^ and A1835 
at ~ 120 Mq yr"^ (see Table 3 in Bildfell et al. and refer- 
ences therein). 

Setting aside the results for MS 1512 for reasons noted 
at the beginning of this Section, Figure [5] indicates that the 
recent star formation in the blue-core BCGs typically has 
an age less than 200 Myrs and contributes mass fractions 
of less than a percent. Using the values of the most likely 
tyc and fyc inferred for the young stellar component in our 
blue-core BCGs, we can estimate an upper limit for the SFR 
in the last ~ 0.2 Gyr for the UV-blue ellipticals by assum- 
ing that a mass of stars fyc x Mt,BCG is assembled over a 
time tyc- The rates span the range ~ 20 — IQOMQ/yr. These 
rates are broadly consistent with the published star forma- 
tion rates determined by other means though we note that 
the agreement is much better for systems with high star for- 
mation rates and that our derived rates tend to be higher 
than the rates derived from spectroscopy for systems with 
low star formation rates. That our results tend to be slightly 
higher is perhaps not a surprise especially when one allows 
for the fact that we have assumed that all of the observed 
UV flux from the blue-core BCGs comes from the stellar 
component whereas this flux may include some (small) con- 
tribution from the AGN. One could also argue that even 
in the absence of any AGN contamination, our estimated 
star formation rates are upper limits because of the way we 
have chosen to model this process as an instantaneous burst, 
as opposed to a prolonged event that extends towards the 
present. While true in general, we do not expect this to be 
an issue here because the age of the young stellar popula- 
tion in our blue-core galaxies is comparable to the typical 
lifetimes of O and B stars. In such cases, the extended and 
instantaneous models should both give similar results. 

4.1 Plausible mechanisms for fueling the recent 
star formation in UV-blue BCGs 

Several of our blue-core BCGs, including A1835 and A2390, 
show evidence for signiflcant accumulation of dust and 
molecular gas in their central regions (Edge et al. 2002, 
Egami et al. 2006). This reservoir is believed to provide the 
fuel for the ongoing star formation but where did this gas 



come from? Three obvious sources are: (1) gaseous flows as- 
sociated with the cooling of the intracluster medium, (2) gas 
deposition during mergers, and (3) gas ejected from the stars 
that comprise the BCG (Menanteau et al. 2001a; 2001b). 

Of these, Bildfell et al (2008) strongly favor the first 
on the strength of the correlations between the presence (or 
absence) of star formation in the BCGs and the global X-ray 
properties for the host clusters. Bildfell et al. (2008) showed 
that the blue-core BCGs are exclusively located in clusters 
that define the high-luminosity edge of the scatter in the 
Lx-Tx plane. In other words, at a fixed temperature, only 
clusters with the highest Lx display observable signatures of 
active star formation. 

This region of the Lx-Tx diagram is known to be pop- 
ulated predominantly by cool-core clusters. In the absence 
of a heating mechanism that can fully compensate for the 
efficient loss of thermal energy through radiative losses, the 
cooling gas would be expected to flow inward towards the 
center of the potential, giving rise to cooling flows. Cool- 
ing flows provide a natural way to account for the reser- 
voir of molecular gas in detected in some of the BCGs. 
More recently, the link between the central cooling timescale 
of the intracluster medium and star formation was further 
strengthened by the results of Cavagnolo et al (2008) and 
Rafferty et al. (2008), who used Chandra observations to 
determine the central gas densities and temperatures for a 
large number of clusters and found that star formation typi- 
cally only occurs in those clusters with central cooling times 
tc < 1 Gyr. Interestingly, a detailed analysis of the individ- 
ual BCG-cluster systems in their CCCP sample led Bildfell 
to discover that while being embedded in a cool core cluster 
is a necessary condition for recent star formation in a BCG, 
it is not a sufficient condition. The BCG must also lie at 
the centre of the cluster potential. It is the combination of 
these two trends that led Bildfell et al. (2008) to support 
the cooling flow hypothesis. 

The fact that all of our blue-core systems show evi- 
dence of a young stellar population (< 200 Myr) suggests 
that either we have caught these systems within 200 Myrs of 
their most recent bout of star formation, or that these sys- 
tems have been forming stars for most of their lifetime. We 
find the latter to be much more plausible especially since all 
the BCGs in the CCCP sample that are hosted by clusters 
with high central X-ray surface brightness and are situated 
within ~ 10h~^ kpc of the center of the host cluster poten- 
tial (identified as the peak of the X-ray surface brightness 
distribution) are blue-core and hence, actively star forming 
systemslfl 

If the star formation in the centrally located BCGs in 
cool core clusters is indeed fueled by cooling fiows, then the 
above-mentioned results also provide insights into the na- 
ture of cooling flows. In the absence of ongoing heating, the 
short central cooling times of observed cool core clusters im- 
ply that radiative cooling would rapidly establish extreme 
centrally-peaked gas density profiles and sharply declining 
central temperature profiles, with gas temperatures tending 
towards zero at the cluster center. However, Chandra and 
XMM-Newton X-ray observations have convincingly demon- 
strated that present-day cool core clusters are not in the 



These systems comprise ~ 20% of the full CCCP systems. 
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Figure 5. Probability maps for tyc and fyc for the UV-blue BCGs. Recall that tyc and fyc represent the age and mass fraction 
of the recent star formation in the galaxy. The recent star formation in the blue BCGs typically has an age less than 0.2 Gyrs and 
contributes mass fractions of less than a percent. 



state just described — their central cores are warmer than 
expected (e.g., Peterson et al. 2003). This has led to claims in 
the literature that the AGNs in the centrally located BCGs 
are injecting sufficient energy into the intracluster medium 
to offset radiative losses and prevent the formation of cooling 
flows. 

We do not dispute that the intracluster gas is being 
heated; rather, based on our results, we assert that once the 
cooling flow is established, it is essentially "on" for most of 
the time and that heating only acts to temper the cooling 
flow, not prevent them. This is contrary to strong feedback 
models of De Lucia et al. (2006). We cannot preclude the 
possibility of temporary disruption by mergers or AGN ac- 
tivity, only that if disrupted, they reform quickly. That is, 
in models where AGN feedback actually manages to heat 
the cluster gas sufficiently so as to halt the flow of cold gas 
towards the cluster center, such a phase must have a rela- 
tively short lifetime (< 200 Myrs), after which the flow and 
the associated star formation resumes. Otherwise, we would 
expect to see a fraction of the centrally located BCGs in 
cool core systems showing "E-|-A"-like characteristics (i.e. 
the latest generation of stars with tyc ^ f Gyr). Admittedly, 
our sample is small and we are mindful of that; however, the 



results for our blue-core systems stand in stark contrast to 
those of the BGCs in non-cool core clusters, in which the 
youngest stellar populations have ages >f Gyr. 

Regarding the other two potential explanations for the 
cold gas and star formation, (Kaviraj et al. 2007a) find the 
recycling of the internal reservoir of gas built up from mass 
loss from the extent stellar population (Menanteau et al. 
200fa; 2001b) is not efficient enough to render the ellip- 
ticals UV-blue. Moreover, if this was the explanation, one 
would also expect the BCGs in non-cool core clusters to 
show evidence of current star formation, at variance with the 
findings of Bildfell et al. (2008) and Rafferty et al. (2008). 
Additionally, one would expect that central star formation 
fueled by recycled gas would tend to preserve and possibly 
even reinforce the radial metallicity gradients that is typi- 
cally found in BCCfS (see Pipino et al. 2006, Pipino et al. 
2008). Cardiel et al. (1998), however, found that in BCGs 
featuring emission-lines (due to recent star formation), the 



A decrease of the metal content in the stars of -0.3 dex per 
decade in radius is typically observed with mild (if any) correla- 
tion with galactic mass in ellipticals (e.g. CaroUo et al. 1993) and 
BCGs (e.g. Fisher et al. 1995, Brough et al. 2007). 
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gradients in the spectral indices are fiat or even positive 
inside the emission-line regions. Outside the emission-line 
regions, and in cooling flow galaxies without emission lines, 
gradients are negative and consistent with those measured 
in giant elliptical galaxies (Fisher et al., 1995). At the same 
time, a residual star formation can lower the central [a/Fe] 
abundance ratio (e.g. Pipino & Matteucci, 2006) and make 
blue BCGs differ from massive ellipticals. Larger samples, 
as the one by Loubser et al. (2008) will shed more light on 
the details of the chemical evolution in BCGs compared to 
normal ellipticals, providing additional clues into the source 
of the gas that is fueling the observed star formation and to 
the integral star formation histories. 

The accretion of gas-rich galaxies by the central BCG 
is another mechanism by which the latter can acquire its 
reservoir of gas. In fact, this mechanism may seem to be 
plausible for BCGs because they are likely to have more 
companions than the average elliptical galaxy in the field 
and because dynamical friction within the cluster halo is 
likely to facilitate the orbital decay of these galaxies. Kavi- 
raj et al. (2007c) find that galaxy mergers of mass ratios 
less than 1:4 and where the cold gas amounts to 20% of the 
merging satellite mass, produce good agreement with the ob- 
served UV colors of regular ellipticals. There are, however, 
several problems with this particular explanation, of which 
the two most important are: (a) The galaxies orbiting in 
or falling through the inner cluster environments are typi- 
cally gas poor; it is believed that they have either lost their 
gas to the intracluster medium via ram pressure stripping, 
or it was consumed by star formation and not replenished 
because of the nature of the local environment, (b) This 
mechanism for delivery of cold gas should, to first order, be 
equally applicable to centrally located BCGs in both cool 
core and non-cool core clusters. Consequently, we would ex- 
pect to see blue-core galaxies in both types of clusters. We 
don't. Additionally, there does not appear to be any out- 
right correlation between presence or absence of a blue-core 
and the local density of galaxies in the neighborhood of the 
BCG. For instance, MS145fJEl has very few satellites but ex- 
hibits a blue core while Abell 2537'' has a red core but many 
satellites. 

Of the three mechanisms discussed above, the cooling 
flow explanation seems to be the most promising primary 
mechanism for fueling the observed star formation. 

Finally, the confirmation of recent star formation in 
early-type galaxies calls for a dedicated survey to detect 
Type II supernovae (SNII) in these systems. From a theoret- 
ical point of view, in fact, the current observational upper 
limit in the central parts of galaxy clusters (0.027 SNu^fl, 
Mannucci et al. 2008) poses a strong constraint on the re- 
cent star formation even if this value is slightly higher than 
for the ellipticals as a whole. In principle, we should be able 
to observe SNII events in blue core BCGs. However, the 
typical brightness of SNII, combined with the fact that the 
recent star formation takes place in the bright centres of 
the galaxies, implies that the probability of finding SNII 
in such galaxies using the current generation of surveys is 
probably low. In fact, according to Mannucci et al. (2007), 



^ Not analysed here, but present in the Bildfell et al sample 
SNuM = SN explosion per century per IO^^Mq 



5-10% of the local SNII remain out of reach of current su- 
pernova surveys. The fraction of missing events rises sharply 
towards z = 1, when about 30% of the SNII might be un- 
detected. Nevertheless, some good candidates do exist, such 
as the BCG in A1835, which have very high star formation 
rates across an extended star forming region, enhancing our 
chances of finding SNII. Furthermore, the SNIa rate seems 
to be slightly higher in radio-loud galaxies, indicating that 
recent star formation may have boosted the prompt channel 
of SNIa progenitors (for details see Delia Valle et al. 2005, 
Mannucci et al. 2006) . 



5 CONCLUSIONS 

In a recent study, Bildfell et al. (2008) showed that the ma- 
jority of Brightest Cluster Galaxies (BCGs) in the Canadian 
Cluster Comparison Project (CCCP) have shallow optical 
colour profiles that become bluer with increasing radius. 
However, a substantial minority of the BCGs in the sam- 
ple (25%) deviates from this simple behaviour, exhibiting 
blue optical cores instead. In this study, we have presented 
a joint analysis of the near ultra-violet (NUV) data from the 
GALEX mission and optical colour profiles for a sample of 
7 BCGs from the Bildfell et al. sample. 

We find that every BCG that has a blue NUV colour 
also shows a blue-core in its optical colour profile. Con- 
versely, BCGs that lack blue cores and show monotonic 
colour gradient are red in the UV. Although the sample 
studied here is small, we demonstrate, for the first time, 
that a one-to-one correspondence between blue cores and 
UV-enhancement is a clear indicator of recent star forma- 
tion. The UV-enhancement is not due to old evolved stellar 
populations such as horizontal branch stars. While we can- 
not outrightly rule out a contribution to the observed UV- 
fiux from an obscured central AGN in the BCG, a number 
of factors including the physical size of the blue cores and 
the fact that obscured AGNs are expected to contribute to 
the optical signal at a few percent level lead us to assert that 
AGNs are not the principal agents of the optical and NUV 
observations. Our assertions are in agreement with other in- 
dependent indicators of star formation in BCGs based on 
optical and infra-red observations (e.g. Cardiel et al. 1998, 
Crawford et al. 1999, Edge 2001, Goto 2005, McNamara et 
al. 2006, Wilman et al. 2006, BildfeU et al. 2008, Cavagnolo 
et al. 2008, O'Dea et al. 2008, Rafferty et al. 2008). 

The young stellar component in all of our blue core 
BCGs typically has an age less than 200 Myrs. The presence 
of ongoing star formation strongly implies that these systems 
have been forming stars for most of their life time, albeit at 
a low rate. Given that all the BCGs in the CCCP sample 
that are hosted by clusters with high central X-ray surface 
brightness and situated within ~ 10h~^ kpc of the center 
of the host cluster potential are blue-core, we must accept 
that all such systems must be steadily forming stars over a 
cosmological time. 

We discuss several possible sources of gas for feeding 
the recent star formation in the blue core BCGs: (a) cool- 
ing flows (Bildfell et al. 2007); (b) recycling of stellar ejecta 
(Menanteau et al. 2001); (c) mergers with gas-rich galaxies 
(Kaviraj et al. 2007c). However, an analysis of each of these 
possibilities combined with unambiguous linkages between 
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the color of the BCG core, the location of the BCG within 
the host cluster, as well as the excess X-ray luminosity (Bild- 
fell et al. 2008) of the cluster and the short central cooling 
time of the intracluster gas (Rafferty et al. 2008) leads us 
to conclude that the cooling flow is the primary mechanism 
for fueling the recent star formation BCGs. In this regard, 
BCGs seem to behave differently from the rest of elliptical 
galaxies. This also means that heating from the AGN largely 
moderates the cooling flow, as opposed to preventing its for- 
mation altogether like in the strong feedback model of De 
Lucia et al. (2006). 

Finally, we note that while we have demonstrated a di- 
rect correspondence between the presence of a blue optical 
core and a blue UV colour in a small sample of BCGs, the 
result needs to be confirmed through an analysis of a much 
large sample of elliptical galaxies, drawn from a wide range 
in luminosity and environments. Among the interesting is- 
sues to consider is whether BCGs in group environments 
because similarly to our cluster BCGs. Future work will fo- 
cus on repeating the analysis presented here on GALEX- 
detected early-type galaxies that have SDSS images (from 
which colour profiles can be extracted) to test the robustness 
of the preliminary results presented in this study. 
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